function r = pulse ( c, w, x )
  r = 1 ./ exp( (x - c) .^ 2 ./ w );
endfunction 
